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Abstract 

The results for the running of the gauge couplings in the MSSM 
are up-dated by proper inclusion of all low scale effects. They are pre- 
sented as predictions for the strong coupling constant in the scenario 
with only two parameters at the GUT scale (ajj and Mjj) and as a 
mismatch of the couplings at the scales ~ 3 x 10 16 GeV and 4 x 10 17 
GeV, when all three couplings are taken as the experimental input. 
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1. The gauge coupling unification QTJ| within the Minimal Supersymmet- 
ric Standard Model (MSSM) has been widely publicized as a successful 
prediction of SUSY-GUTs [2-6]. It is also often discussed in the context of 
stringy unification, with M ST ~ 4 x 10 17 GeV f7|. In this paper we up-date 
the results for the running of the gauge couplings in the MSSM by proper 
inclusion of all low energy effects such as the best precision of the input pa- 
rameters at the electroweak scale and the non-logarithmic contribution from 
the superpartner thresholds R |9|, [10, |TTJ . 

The unification idea is predictive with respect to the behaviour of the 
SU(3) x SU(2) x U{1) gauge couplings if physics at the GUT scale can 
be described in terms of only two parameters: ot\j and M\j (minimal 
unification). Then we can predict e.g. a s (M z ) in terms of (^em{Mz) and 
siia 2 9 w(M z ) (it is worth remembering that sin 2 9 W (M Z ) and a s (M z ) are 
at present known with 0.1% and 10% accuracy, respectively). More precisely, 
the prediction for the strong coupling constant in addition depends on the 
superpartner spectrum which will, hopefully, be known from experiment. For 
now, these are free parameters and, denoting them globally by Tsusy (see 
the discussion in section 2) we get 

a s {M z ) = F(sm 2 9 w (M z ),a EM (M z ), T SUSY ) (1) 

This approach may, however, be too restrictive as it is generally expected 
that there are non- negligible GUT/string threshold corrections to the run- 
ning of the couplings (such as heavy threshold and higher dimension oper- 
ator effects). Then, strictly speaking, all predictivity is lost. However, it is 
still very interesting to reverse the problem: take the values of all the three 
couplings at M z as input and use the bottom-up approach to study the 
convergence of the couplings in the framework of the MSSM. With the same 
precision calculation and as a function of the SUSY spectrum one can, then, 
discuss the mismatch of the couplings at any scale of interest and for any 
value of a s (M z ), within its 10% experimental uncertainty. It is convenient 
to introduce the "mismatch" parameters at scale Q: 

A(Q) = a ' W) 7 ° 2(Q) (2) 

and 

MQ) = - — TTyT (3) 

(the latter are directly related to large scale threshold corrections). Of par- 
ticular interest are D^(Mu), where My is defined as the scale of unification 
of the SU(2) x U(l) couplings (i.e. the scale at which D% = A\ = 0), and 
Di(Msr), i = 1,3, with M ST = 4 x 10 17 GeV 0, 0] (and corresponding 
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Ajs). Clearly, we get this way constraints on physics at the high scale, if 
it is supposed to have unification and the MSSM as the low energy effective 
theory. We can also read this information as a hint whether the latter two 
assumptions look plausible. 

In this paper we present our results both as the prediction for a s (M z ) 
in the minimal unification scenario and as a prediction for the mismatch 
parameters at M v and M ST , as a function of a s . 



2. We begin with the discussion of the experimental information. Let 
us first suppose that the (non-supersymmetric) SM is the correct effective 
theory at the electroweak scale. In this theory the couplings g 3 , g 2 , g\ of 
the SU(3) x SU(2) x SU(1) gauge groups , at M z and in the MS scheme 
are usually quoted as the values of a em 

(M z ), sin 2 6 W {M Z ) and a s {M z ). 
The electromagnetic coupling constant and the Weinberg angle in the SM 
are now known with very high precision. The value of oiem{Mz) * s obtained 
from the on-shell ot% M = 

1/137.0359895(61) via the 1-loop RG improved 



relation [13 



a EM (M z ) = ( 4 ) 
1 — A« 



where 



Aa = 0.0682 ± 0.0007 +— log — ^ log — - 5 

2tt 6 M z 9vr 6 M z v ' 

The main uncertainty comes from the continuous hadronic contribution to 
the photon propagator. We explicitly show the top quark mass dependence 

of a EM (M z ). 

The most precise value of sin 2 9w(M z ) in the MS scheme is at present 
obtained in terms of Gf, M z and cxem- The result depends on m t and 
M^o (the top and the SM Higgs boson masses respectively) and to a very 
high precision is given by the following effective formula [] : 

sin 2 9 W (M Z ) = 0.23166 ±0.0003 + 5.4 x 10~ 6 h - 2.4 x 10~ 8 /i 2 

- 3.03 x 10" 5 t - 8.4 x 10~ 8 t 2 (6) 

where h = M^o — 100 and t = m t — 165 (both masses in GeV). The main 
source of the error is again the hadronic uncertainty in the photon propagator. 
E.g. for m t = 180 GeV and M^o = 100 GeV we get sin 2 #iy(M z ) = 0.2312. 

The value of a s (M z ) is known with much worse precision and depending 
on the method of determination, the values in the range 0.11-0.13 are quoted 

1 Notice the small change as compared to ref. || which is due to the inclusion of QCD 
corrections according to ref. |0. Another, frequently used, fit is given in ref. j| and its 
update in H. 
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It is interesting that the lower part of this range is favoured by low 



energy determinations of a s ||T6| , |17| and by a fit to all electroweak data in 



the framework of the MSSM 19 



Once (yfjS at Mz in the SM are extracted from the data, the 2-loop RGE 
can be used to get them at higher scales. Passing through the thresholds of 
superpartners the running of the couplings is subject to subsequent modifi- 
cations of the /3-functions || with, finally, MSSM RG equations above all 
the thresholds. Treating the threshold corrections at the 1-loop level (con- 
sistently with the 2-loop RGE) this procedure gives 0: 



1 1 C * 9^TAh Mk 



+ 2£A^log^ (7) 

k % 



+ 2bf SSM log — + two - loop contribution 

M 7 



where we have made explicit the MS — > DR conversion factor with C\ = 
0, C2 = 2, C3 = 3 [pO[l . Mk are the superpartner masses and Abu- are 
their contributions to the one-loop j3 functions of the couplings ctj. This 
is the correct result for the running of the gauge couplings at the two-loop 
accuracy as long as the contribution to the SM from the (non-renormalizable) 
higher dimension operators, left over after decoupling of superpartners, can 
be neglected in the process of extracting g^Mz) from the data (we shall call 
it the Leading Logarithmic Threshold (LLT) approximation). This requires 
M k 3> M z for all superpartner masses. 

Assuming that there are only two GUT scale parameters: au and 
Mu (i.e. assuming that the potential GUT scale corrections to the gauge 
coupling unification are negligible) we can predict one of the couplings at 
Mz scale, e.g. a s [Mz) in terms of the other two and of the superpartner 
masses, which are at present free parameters. In the LLT approximation the 
dependence of the prediction for a s (Mz) on the supersymmetric spectrum 
can be described by a single effective parameter T S usy @ : 

a s {M z ) = f(a 1 (M z ),a 2 (M z ),T SUSY ) (8) 
= f(G F ,M z ,aEM( y M z ),m t ,M (t> o,TsusY) 



where 111811 : 



j_ 

V.) 



2 af M 



{Mz) differs from a s {Mz) by a threshold correction from the top quark: a 3 1 



aj x + ±log(m t /M z ). 
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The effective parametrization in terms of Tsusy is exact for one-loop RGE 
and the correction due to the superpartner spectrum then reads: 



19 , Tsusy 



af M ~ a° ' 2vr 14 g M» 



(a® is the value predicted without the inclusion of threshold corrections). 
With two-loop equations there is some (weak) dependence on the details 
of the spectrum through the dependence on the spectrum of the two-loop 
contribution on the way up to M v . 

The prediction of the eqns. (|7HlO|) may be subject to important correc- 
tions if some of the superpartner masses are 0(M Z ). Then the renormaliz- 
able SM is not the correct effective theory at the electroweak scale and the 
non-renormalizable terms should be included when extracting the couplings 
from the data. Equivalently, we can work at M Z in the framework of the 
full MSSM, extract from the data the MSSM couplings including full 1-loop 
threshold contribution from SUSY loops (not just the leading logarithms) 
and study the unification of the MSSM couplings. (Note that in the LLT 
approximation an equivalent interpretation of equation (H) is: 



: = : + 2b MSSM be — ( 1 1 ) 

+ two — loop contribution 



with 



af SSM {M z ) ~ a? M (M z ) 12vr + 2 ^ Akk hg M z (12) 

where RG running with the MSSM f3— functions starts directly from M z , 
and the threshold corrections are absorbed in a redefinition of gi(M z ) SM — > 
g t (Mz) MSSM ) 

The outlined program has been accomplished by several groups: ||, |10 



[DJ . Clearly, the values of the MSSM couplings extracted from the data 
depend now on the superpartner masses Mk, e.g.: 

sin 2 9 W (M Z ) MSSM = f(Gp, M Zl a EM , m u M k , M^o , tan (5) 

not only by logarithmic terms as in eq.(|l"2"D, but also by terms 0(M Z / M^) 
and the additional corrections may be ~ 1% for sin 2 9 W (M Z ) MSSM as 
shown in ref. [[|. 

In our analysis we also use a^ SSM with the oblique non-logarithmic 
corrections included P, H but they are unimportant for generic spectra which 



3 Also Map must be replaced by M^o and tan j3. 
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have coloured sparticles rather heavy. We also use the properly extracted 



a 



MSSM 
EM 



The impact of the non-leading SUSY corrections on the prediction for 
a s (Mz) is illustrated in Fig. 1 for a generic sparticle spectrum obtained in 
the minimal supergravity model (with universal boundary conditions for the 
soft SUSY breaking scalar mass parameters at the GUT scale) with radiative 
electroweak breaking and squark masses below 2 TeV [|2T|. The results for 
a s (M z ) are plotted as a function of T SUS y defined in eq. (Q). We 
compare the results obtained in the LLT approximation for the superpartner 
thresholds, eq. (0), with their complete inclusion at the one-loop level, as in 
ref. 0. The non-leading corrections increase the predicted value of a s (M z ) 
for T S usy < 100 GeV. We conclude that unification without GUT threshold 
corrections predicts e.g. for m t = 160 (180) GeV a s {M z ) > 0.126 (0.128) 
for T SUSY < M z and a s {M z ) > 0.121 (0.123) for T SUSY < 300 GeV. 
It is clear from eq.(^) that T S usy depends strongly on /i, m^, m g and 
weakly on the other SUSY masses. In models with the GUT relation 

M 3 (M V ) =M 2 (M V ) , (13) 

to a very good approximation: 

W ~H^)J ~^ (14) 

and large Tsusy means very large higgsino mass. From the naturalness of 
the Higgs potential ^2], |j it follows then that also the other sparticle masses 



are to be heavy. For instance in the generic spectrum obtained in models 
with radiative breaking and universal boundary conditions for the soft scalar 
masses at the GUT scale T S usy = 300 GeV corresponds to the squark 
masses 0(2 TeV). Of course, large values of Tsusy can be obtained 



also for small /i with a spectrum which violates the GUT relation (0), i.e. 
with a large ratio m^/m g . However, it is very difficult to imagine such a 
scenario without losing the motivation for the minimal unification itself. 

3. The assumption about negligible GUT scale corrections to coupling 
unification may be too restrictive. Various groups have discussed the GUT 
threshold corrections (dependent on the GUT model) ]24|, |25|, [|, [| and 



0(Mu/Mpi) corrections |27| , |28|j . Admitting non-negligible but strongly 
model dependent GUT scale corrections means that, strictly speaking, the 
predictivity is lost and one can only use the bottom-up approach: measure 
gi(M z ) with better and better precision, measure the sparticle spectrum and 
study the convergence of the couplings at the large scale. Useful mismatch 

4 The uncertainty in a s induced by the errors in eqs. (^,||) is ~ 0.0015 Q. 
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parameters then are: D 3 (Mu) and A 3 (Mu) (eqs. (0,0); Mu is defined by 
unification of ai and oc^). A 3 (Mu) is directly related to the GUT threshold 
corrections. Again, neglecting O(MzfMk) non-renormalizable terms, both 
D 3 (Mu) and A 3 (Mu) are functions of a i: i = 1,2,3 and the effective 
Tsusy- Inclusion of non- leading supersymmetric threshold corrections 
brings in additional dependence on the spectrum with, however, T S usy still 
a useful parameter to present the results. We show them as a function 
of Tsusy i n the LLT approximation for the SUSY thresholds and with 
their complete inclusion in Fig. 2 for our generic spectra for m t = 180 GeV , 
tan/3 = 10 and for three values of a s =0.11, 0.12, 0.13. In Fig. 3 we plot 
the same mismatch parameters as a function of a s for our generic spectra. 
The LLT results for two fixed values of T S usy = 300 GeV and 1 TeV are 
also shown for comparison. 

The general conclusion is that in the range a s (Mz) =0.11 — 0.13 and 
Tsusy = (20 — 10 3 ) GeV the gauge couplings do unify within the accuracy 
better than 7% for m t = 160 GeV and 0(8%) for m t = 180 GeV, with the 
maximal mismatch for low values of a s (Mz) and Tsusy- Is this mismatch 
a lot or a little depends on the GUT model and the expected magnitude of 
the GUT scale corrections in it 



4. In stringy unification the unification scale is no longer a free parameter. 
It is related to the value of the unified coupling JT3]: 

M ST = 9st x 5.27 x 10 17 GeV ~ 0(4 x 10 17 GeV) (15) 

It is interesting to study within the bottom-up approach the mismatch pa- 
rameters D 3 , D 1} A 3 , Ai (eqs. at the scale M ST = 4 x 10 17 
GeV. The results as a function of a s {Mz) are shown in Fig.4. We use 
again our sample of generic spectra. The results for very heavy spectra with 
Tsusy = 1 TeV and 5 TeV obtained within the LLT approximation are 
also shown. The general conclusion which can be drawn from these plots is 
that the mismatch of the couplings a 3 and a 2 as well as a± and a 2 at 
Mst is > O(10%). Therefore, to achieve unification, the string threshold 
corrections have to be large at the string scale and in addition must conspire 
so that they are small at the GUT scale, i.e. that the approximate unifica- 
tion occurs at Mu ~ 3 x 10 16 GeV. It is also worth pointing out that the 
dependence of Di(Mst) and D 3 (Mst) on the supersymmetric spectrum 
is different and the spectrum which diminishes the first enhances the second. 

It is possible to take the attitude that the value of ai at the string 
scale is unconstrained because the Kac-Moody level of the U(l) group can 
be treated as a free parameter k\. [] We then have: 

kioti = a 2 = a 3 (at M S t) (16) 

In our convention k% = 1 in the case of SU(5) - type unification. This differs from 
the definition adopted in ref. 0. 



24, 25, Ea 11]. 
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In this case (and for negligible stringy threshold corrections) our parameter 
D\ is related to the parameter k\\ 



ki = • (17) 

With our generic spectra we get: k\ = 0.88 — 0.92. However, even then 
we are still faced with a large mismatch between «3 and «2 which for our 
generic spectra requires large string threshold corrections. 

Finally it is interesting to go beyond the discussion based on our generic 
spectra and to address the following two questions: 

1) Does there exist a pattern of the MSSM spectrum which shifts the 
unification point of all three couplings to M$t with negligible stringy thresh- 
olds? 

2) Suppose ki 7^ 5/3 (i.e. a>\ 7^ ct2 at M$t) and helps to unify 
cti and a.2- Are there MSSM spectra which unify and ai at M$t with 
negligible stringy thresholds? 

In order to answer these questions it is useful to introduce two new effec- 
tive parameters describing the impact of the SUSY spectrum on unification 
of «i and «2 and a% and «3 separately. From eqs. (|TTD and (|12"D we 
have: 



M V = M°(^-X (18) 

X 1 SUSY J 

where M v (M£) is the crossing point of a\ and with SUSY threshold 
corrections included (neglected) and 



T'susr = (M*o A$) 25 / \ , (19) 



All generations have the same masses in the above formula but a general- 
ization is straightforward. Noticing that all (none) of the sparticles in the 
denominator (numerator) are SU(2) singlets one can write a simplified for- 
mula: 

M 2 

T ' SUSY = Mb (20) 

with obvious definitions of the averages Ml and Mr. In Fig. 5 we plot Mu 
resulting from formula ([18]). For M v = 4 x 10 17 we need T' SUSY = 2 x 10 -2 
GeV which means that for Ml = Mz we would have Mr = 400 TeV. As 
we can see the answer to question 1 is negative. Regardless of the a 2 — a 3 
unification, bringing Mu up to Mst would require an unacceptable Mr. 
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Turning to question 2 we study the correction to af (Mz) predicted 



from the condition « 3 (Mst) = 0:2 (Mst) 
From (|II|) and 0) we obtain: 



induced by the SUSY thresholds. 



1 



a 



SM 



1 1 19, Md 
-« + 2-^ l0g MS 



Where 



M 



D 



M s 



M 



A 



M, 



a 



(21) 



m 



(22) 



and 



a.. 



is the value predicted without the inclusion of SUSY threshold 



corrections. In Fig. 6 we show predicted with the use of formula (pl|) 
(with a° obtained from the two-loop RGE) for M ST = (3.5 , 4.0 , 4.5) xlO 17 
GeV as a function of the ratio Mr, /Ms- In order to get ot s < 0.13 one 
needs Md/Ms > 20. At this point we disagree with the recent analysis of 

a s (M z ) 



ref. [ 2S 



which reconciles a s (Mz) = 0.118 with the stringy unification for 
the SUSY spectra with smaller hierarchies. Taking masses in the numerator 
of our eq. (|22| ) ~ 30 times larger than masses in the denominator we still 
get a s > 0.125 as can be seen from Fig. 6 (]. This disagreement is mainly 
due to the use of one-loop RGE in ref. 
2 9w(Mz) (more appropriate for 



29| and somewhat higher value of 



m 



ji/m 



m t ~ 160 GeV). The ratio M D /M S is 
which in the case of string unification is 



sin owymz 

dominated by the ratio 

more model dependent than for GUTs |3(], [H| • I n particular it is conceivable 
in this case that > rrig. For Ms = Mz we would get Md > 2 

TeV but in fact for so light spectrum the non-logarithmic corrections could 
raise the predicted value of a s as is evident from Fig.l. For Ms = 150 
GeV, when non-logarithmic effects are small, we get Md > 3 TeV. We 
conclude that it is possible to raise the a 2 — a 3 unification scale up to M S t 
but only with highly unnatural SUSY spectra, with the heaviest sparticles 
above 3 TeV and with a s (Mz) — 0.13. Otherwise large string threshold 
corrections are needed. 



5. We have discussed the impact of SUSY thresholds on the unification 
of gauge couplings in the framework of GUT and string theories. Non- 
logarithmic SUSY corrections can be important for the phenomenologically 
interesting case of light superpartners. These corrections always reduce the 
value of siia 2 9w(Mz) MSSM which in turn raises the value of a s (Mz) SAI 
predicted from SUSY unification. In the minimal unification scenario (i.e. 
with negligible GUT scale corrections to the running of the couplings) one 
gets a s (M z ) > 0.121(123) for the effective parameter T SUS y < 300 GeV 

6 For a s (M z ) — 0.12 the masses should be split by a factor of at least 60. 
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and a s (M z ) > 0.115(117) for T SU sy < 1 TeV, for m t =160 and 180 
GeV, respectively. For the generic spectra in the minimal supergravity model 
Tsusy ~ 1 TeV corresponds to very heavy sfermions, e.g. squark masses are 
(9(5 TeV). More generally in the bottom-up running the couplings do unify 
within a few percent accuracy even for low a s [Mz) and small values of 
Tsusy, e.g. the mismatch between a 3 and a 2 at the scale of unification 
of ai and a 2 is generically below (9(5%). The mismatch of the couplings 
at M S t = 4 x 10 17 GeV is much larger, typically (9(10%) or more, and it 
cannot be eliminated by any sensible superpartner spectrum. String unifica- 
tion requires, therefore, large string threshold corrections (which, however, 
may not be unrealistic |3!|) which conspire to give the effective unification 
scale ~ 3 x 10 16 GeV. The scenario with «i and a 2 unified by treating 
the Kac-Moody level k\ as a free parameter is not particularly helpful with 
regard to the coupling unification at M$t (and is rather uneconomical). 
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FIGURE CAPTIONS 



Figure 1. 

a s (Mz) predicted by the minimal unification as a function of Tsusy for 
different values of m t and tan (3. Squares (stars) correspond to the LLT 
(full) calculation of the supersymmetric thresholds for a generic sample of 
SUSY spectra obtained in the minimal supergravity model with radiative 
electroweak breaking and universal boundary conditions. 

Figure 2. 

Mismatch parameters A3 (My) and D 3 (Mjj) for m t = 180 and tan/5 = 10 
as a function of T S usy for several values of a s . Squares (stars) correspond 
to the LLT (full) calculation of the supersymmetric thresholds. Sample of 
spectra as in Fig.lc. Solid lines extrapolate A3 (My) and D$(Mu) in the 
LLT approximation up to Tsusy = 1 TeV. 

Figure 3. Mismatch parameters A 3 (Mj/) and D 3 (Mjj) as a function of 
a s for the same sample of spectra as in Fig.lc. Squares, stars and circles 
show the results of the full calculation for spectra with Mq < 500 GeV, 
500 GeV < M Q < 1 TeV and 1 TeV < Mq < 2 TeV. respectively. For 
comparison the LLT calculation for a spectra with Tsusy =300 GeV (1 
TeV) are marked by the solid (dashed) lines. 

Figure 4. 

Mismatch parameters A3, D%, A% and D% at the string scale M$t = 
4 x 10 17 as a function of a s for our generic spectra. Markers as in Fig.3. 
Solid, dashed and dash-dotted lines correspond to the LLT calculation for 
Tsusy =300 GeV, 1 and 5 TeV respectively. 

Figure 5. Mu as a function of T' SUSY from eq. (jl8|) for m t = 180 GeV. 
Figure 6. 

Prediction for a s from the condition a^Msr) = ^(Mst) with Mst = 3.5, 
(solid) 4.0 (dashed) and 4.5 x 10 17 GeV (dash-dotted) with the use of 
eq. (PH) as a function of the parameter Mm/ Ms for m t = 180 GeV. 
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